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[i]  A  new  mean  dynamic  topography  (MDT)  for  the  Bering  Sea  is  presented.  The  product 
is  obtained  by  combining  historical  oceanographic  and  atmospheric  observations  with 
high-resolution  model  dynamics  in  the  framework  of  a  variational  technique.  Eighty 
percent  of  the  ocean  data  underlying  the  MDT  were  obtained  during  the  last  25  years  and 
include  hydrographic  profiles,  surface  drifter  trajectories,  and  in  situ  velocity  observations 
that  were  combined  with  National  Centers  for  Environmental  Prediction  (NCEP)/National 
Center  for  Atmospheric  Research  (NCAR)  atmospheric  climatology.  The  new  MDT 
quantifies  surface  geostrophic  circulation  in  the  Bering  Sea  with  a  formal  accuracy  of 
2-4  cm/s.  The  corresponding  sea  surface  height  (SSH)  errors  are  estimated  by  inverting 
the  Hessian  matrix  in  the  subspace  spanned  by  the  leading  modes  of  SSH  variability 
observed  from  satellites.  Comparison  with  similar  products  based  on  in  situ  observations, 
satellite  gravity,  and  altimetry  shows  that  the  new  MDT  is  in  better  agreement  with 
independent  velocity  observations  by  Argo  drifters  and  moorings.  Assimilation  of  the 
satellite  altimetry  data  referenced  to  the  new  MDT  allows  better  reconstruction  of 
regional  circulations  in  the  Bering  Sea.  Comparisons  also  indicate  that  MDT  estimates 
derived  from  the  latest  Gravity  Recovery  and  Climate  Experiment  geoid  model  have  more  in 
common  with  the  presented  sea  surface  topography  than  with  the  MDTs  based  on  earlier 
versions  of  the  geoid.  The  presented  MDT  will  increase  the  accuracy  of  calculations  of  the 
satellite  altimeter  absolute  heights  and  geostrophic  surface  currents  and  may  also  contribute 
to  improving  the  precision  in  estimating  the  geotd  in  the  Bering  Sea. 

Citation:  Panteleev,  G.,  M.  Yaremchuk,  P.  J.  Stabeno,  V.  Luchin,  D.  A.  Nechaev,  and  T.  Kikuchi  (201 1),  Dynamic  topography 
of  the  Bering  Sea,  J.  Geophys .  Res.,  116 ,  C05017,  doi:  10.1 029/20 10JC006354. 


1.  Introduction 

[:]  The  density  and  diversity  of  oceanographic  observa¬ 
tions  in  the  Bering  Sea  increased  significantly  during  the  last 
decades.  The  observations  comprise  conventional  tempera- 
turc/salinity  data,  a  large  number  of  high-quality  velocity 
time  series  from  moorings  [e.g.,  Woodgate  et  al .,  2005],  and 
surface  [ Niiler ,  2001]  and  subsurface  (Argo)  profiling  floats 
[Wilson,  2000].  The  surface  drifter  program  has  allowed  a 
quantitative  determination  of  the  basic  features  of  the  Bering 
Sea  surface  circulation  [Stabeno  et  al .,  1999]. 

[3]  Since  1992,  satellite  radar  altimetry  has  become  a 
conventional  tool  for  remote  monitoring  of  global  sea  level 
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variations.  Development  of  instrumental  technology  and 
processing  techniques  and  utilization  of  multisatcllitc  data 
sets  have  reduced  satellite  altimetry  product  errors  to  the 
point  that  these  products  can  be  used  to  detect  sea  surface 
height  (SSH)  variations  associated  with  ocean  currents  and 
to  resolve  upper  ocean  mesoscale  eddies.  Combined 
together,  data  from  the  TOPEX/Poscidon,  Jason- 1,  Jason-2, 
Earth  Resources  Satellite  (ERS)-l  and  -2,  Envisat,  and 
Geosat  Follow-On  (GFO)  missions  span  the  time  period 
from  October  1 992  to  present. 

[4]  In  contrast  to  the  increasing  accuracy  of  detecting  SSH 
variations,  the  difference  between  a  time-averaged  sea  sur¬ 
face  and  geoid  (mean  dynamic  topography  (MDT))  can  be 
retrieved  from  altimetry  with  much  lower  precision.  The 
problem  is  due  to  the  large  uncertainties  in  the  geoid  models. 
Although  the  recent  Gravity  Recovery  and  Climate  Experi¬ 
ment  (GRACE)  mission  (http://www.csr.utexas.cdu/grace/) 
significantly  improved  the  geoid,  it  remains  too  coarse  to  be 
directly  used  for  circulation  studies  in  the  marginal  seas.  The 
situation  may  improve  w  ith  the  recent  launch  of  the  Gravity 
field  and  steady  state  Ocean  Circulation  Explorer  (GOCE) 
mission. 

[5]  In  the  meantime,  a  number  of  research  groups  have 
developed  methods  to  combine  various  data  with  altimetry 
to  obtain  more  accurate  estimates  of  the  global  MDT  [Niiler 
et  a l .,  2003;  Rio  et  al .,  2005;  Maximenko  et  al. ,  2009].  For 
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various  reasons,  all  these  products  have  substantial  defi¬ 
ciencies  in  the  Bering  Sea.  The  common  cause  is  that  sta¬ 
tistical  assumptions  underlying  these  global  products  were 
not  fine-tuned  locally,  but  SSH  statistics  strongly  depends 
on  local  dynamics,  which  is  affected  in  turn  by  stratification 
and  smaller-scale  topographic  features. 

[6]  Alternative  approaches  include  MDT  estimation  as  an 
ensemble  average  of  regional  model  solutions  [c.g., 
Bingham  and  Haines ,  2006],  or  as  a  diagnostic  solution  of 
the  ocean  model  forced  by  seasonal  climatology  [Foreman 
et  al.,  2008].  Both  methods  take  into  account  dynamical 
information  but  have  certain  limitations  in  assessing  the 
MDT  errors  due  to  limited  ensemble  size  and  uncertainties 
in  model  forcing. 

[7]  The  four-dimensional  variational  (4DVar)  data 
assimilation  technique  computes  MDT  as  a  component  of 
the  mean  climatological  circulation;  this  method  satisfies 
model  equations  on  the  one  hand,  while  on  the  other  hand  it 
provides  the  least  discrepancy  with  observations.  The 
4DVar  technique  has  proved  to  be  a  useful  and  efficient  tool 
in  numerous  ocean  circulation  studies  [e.g.,  IVunsch,  1996; 
Awaji  et  al .,  2003;  Panteleev  et  al.,  2006a].  A  disadvantage 
of  this  approach  is  its  computational  cost  which  prevents 
production  of  global  4DVar  analyses  [Stammer  et  al .,  2002; 
Menemelis  et  al. ,  2009]  at  resolutions  better  than  0.25°. 

[s]  On  the  basin  scale,  variational  inversions  at  high  lati¬ 
tudes  are  known  in  literature  for  almost  two  decades.  The 
early  efforts  were  limited  to  3Dvar  steady  state  problems  at 
modest  resolutions  with  simplified  dynamics  [ Brasseur , 
1991;  Grotov  et  al .,  1998;  Yaremchuk ,  2001].  More  recent 
four-dimensional  variational  inversions  were  made  at  0.1°- 
0.2°  resolution  with  primitive  equation  constraints  [e.g., 
Panteleev  et  al .,  2006a,  2010;  Powell  et  al .,  2008;  Hoteit 
et  al .,  2010].  Panteleev  et  al.  [2006b]  reconstructed  the 
mean  summer  climatological  circulation  in  the  northern 
Bering  Sea  and  found  a  reasonably  good  agreement  with  the 
Argo  drifters  velocities  at  1000  m  [Yoshinari  et  al.,  2006]. 

[y]  In  the  present  study  we  employ  a  4DVar  inversion  of  a 
primitive  equation  model  to  reconstruct  the  MDT  in  the 
Bering  Sea.  The  model  has  a  horizontal  resolution  of  18  km, 
which  is  slightly  larger  than  the  average  baroclinic  radius  of 
deformation  (15  km  [c.g.,  Chelton  et  al .,  1998]).  This  choice 
both  insures  stability  of  the  tangent  linear  and  adjoint 
models  and  allows  the  most  important  topographic  features  to 
be  resolved,  including  the  passes  through  the  Aleutian  Arc. 

[10]  Special  attention  is  paid  to  a  posteriori  error  analysis 
of  the  optimal  SSH  field,  which  is  accomplished  by  in¬ 
verting  the  Hessian  matrix  of  the  assimilation  problem.  This 
inversion  is  performed  in  the  low -dimensional  subspace 
spanned  by  the  leading  modes  of  the  SSH  variability 
observed  by  satellites.  To  validate  the  product,  we  also 
estimate  its  properties  in  selected  regions  of  the  Bering  Sea 
against  the  available  data  and  compare  it  with  the  global  high- 
rcsolution  MDTs  derived  recently  from  in  situ  observations, 
altimeter  data,  and  geoid  models.  These  comparisons  indi¬ 
cate,  in  particular,  that  the  latest  GRACE  geoid  model  pro¬ 
vides  an  improved  MDT  estimate  in  the  open  sea  regions. 

[11]  The  paper  is  organized  as  follows.  In  section  2  we 
describe  the  data  sets  used  and  the  techniques  of  data 
assimilation  and  error  analysis.  The  optimized  MDT  and 
results  of  climatological  reconstruction  of  the  Bering  Sea 
circulation  are  presented  in  section  3.  Also  in  section  3  local 


circulations  in  the  Amukta  Pass,  Eastern  Bering  Sea  shelf, 
the  Bering  Slope,  and  the  Kamchatka  and  Alaska  currents 
are  reconstructed  using  four  different  MDTs  (including  the 
present  one)  and  the  respective  flow  fields  are  compared 
against  independent  drifter  and  mooring  data.  In  section  4 
the  results  are  summarized. 

2.  Methodology 

2.1.  Data 

[12]  The  following  data  were  utilized  in  the  reconstruction 
of  the  MDT: 

[13]  1.  The  81,91 1  temperature/salinity  profiles  that  were 
collected  in  the  Chukchi  and  Bering  seas  between  1932  and 
2004  (Figure  1 ).  This  database  includes  bottle  data,  mechan¬ 
ical  bathythermograph  data,  high-resolution  conductivity- 
temperature-depth  (CTD)  profiles,  expendable  bathythermo¬ 
graph  and  PALACE  Argo  float  data.  More  than  50%  of  these 
data  were  obtained  during  the  period  1980-2004.  The  major 
part  of  the  data  was  obtained  from  the  Russian  archives  in 
R1HM1-WDC  (http://nodc.meteo.ru/nodc/),  JODC  (http:// 
www.jodc.go.jp/).  University  of  Alaska  (http://www.ims.uaf. 
edu).  databases  of  the  World  Ocean  Data  Center  [Conkright 
et  al.,  2002],  and  the  Argo  Global  Data  Assembly  Centre 
(http.V/www.coriol  is.cu.org/).  Preprocessing  of  the  data 
included  quality  control,  which  consisted  of  averaging  the 
profiles  in  time  and  over  the  28  *  28  km  bins  and  computing 
the  standard  deviations  crT ,  crs  from  statistics  w  ithin  the  bins. 
Spatial  distributions  of  the  resulting  mean  temperature  T 
(Figure  2a),  salinity  S ,  aT ,  and  as  were  used  in  the  recon¬ 
struction  of  the  mean  Bering  Sea  state.  The  values  of  <77-  and 
as  varied  within  the  ranges  1 .5°C^4°C  (Figure  2b)  and  0.1- 
2.0  ppt  near  the  surface  and  decreased  to  0.1  °C  and  0.03  ppt 
respectively,  in  the  deeper  layers.  The  total  number  of 
hydrographic  data  points  used  in  the  analysis  was  184,109 
for  temperature  and  178,529  for  salinity. 

[u]  2.  About  500  satellite-tracked  drifter  trajectories  from 
the  Fisheries  Oceanography  Coordinated  Investigations 
(FOCI)  database  (http://ww,w.pmcl. noaa.gov/foci)  and  84 
drifter  trajectories  from  the  Global  Drifter  Program  (GDP) 
database  (http://wrWAV.aoml. noaa.gov/phod4iac0  observed 
during  1984-2004.  The  FOCI  surface  drifters  were  drogued 
at  40  m  and  GDP  drifters  had  drogues  at  15  m.  Preliminary 
analysis  of  these  data  included:  (1)  temporal  low-pass  filter¬ 
ing  of  the  trajectories  with  7  day  cutoff  period,  and  (2)  cal¬ 
culating  the  mean  gridded  velocities  and  the  corresponding 
error  variances  through  spatial  averaging  of  the  drifter 
velocities  within  the  30  km  circles.  The  average  error  vari¬ 
ance  ranged  from  5  to  20  cm/s.  Only  55 1 5  gridded  velocities 
obtained  from  averaging  of  at  least  three  different  surface 
drifters  were  assimilated  (Figure  3a).  Most  (82%)  of  these 
velocity  data  were  obtained  from  FOCI  drifters  driven  by 
currents  below  the  Ekman  layer. 

[15]  3.  Velocity  time  series  from  57  moorings.  Most  of 
these  data  come  from  the  Alaska  Ocean  Observing  System 
(AOOS)  database  (http://www.aoos.org)  encompassing  the 
period  of  1970-2005.  Many  of  the  AOOS  velocity  time 
series  are  only  2-4  months  long  and  cannot  be  a  source  for 
reliable  estimates  of  climatological  currents.  To  avoid  con¬ 
tamination  of  the  optimal  solution  by  a  presumably  strong 
“subseasonaP  signal,  error  variances  of  the  AOOS  monthly 
mean  data  were  taken  as  the  largest  of  the  following  three 
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Figure  1.  Coverage  of  the  model  domain  (solid  rectangle)  by  hydrographic  stations.  Bathymetry  con¬ 
tours  of  1000  and  3000  m  are  shown. 


values:  5  cm/s  (an  estimate  of  seasonal  variability  of  the 
vector- averaged  velocity  amplitude),  or  20%  of  the  monthly 
mean  velocity  amplitude,  or  the  RMS  variation  of  the 
original  velocity  time  series.  The  total  number  of  assimilated 
velocity  observations  from  this  source  was  1 14  (Figure  3b). 

[16]  4.  The  Bering  Strait  transport  estimate  of  0.9  ±  0.2  Sv 
was  taken  from  Woodgate  et  al  [2005].  The  estimates  of  this 
transport  arc  based  on  the  data  collected  since  1990-2004. 

[17]  5.  The  wind  stress  and  surface  heat/salt  fluxes  were 
taken  from  the  National  Centers  for  Environmental  Predic¬ 
tion  (NCEP)ZNational  Center  for  Atmospheric  Research 
(NCAR)  climatologies  averaged  over  the  period  of  1948- 
2006  (http://www.csrl.noaa.gov/psd/data/griddcd/data.ncep. 
reanalysis.derived.html).  Since  these  products  may  contain 
substantial  errors  in  the  Bering  Sea  [Ladd  and  Bond ,  2002], 
we  used  wind  stress  and  surface  flux  data  with  relatively 
high  error  variances  equal  to  40%  of  their  spatiotcmporal 
RMS  variation  over  the  basin.  The  total  number  of  surface 
flux  observations  was  34,888. 

[is]  All  the  above  data  sets  differ  in  temporal  coverage;  in 
addition,  their  temporal  distribution  is  uneven.  We  estimated 
that  80%  of  all  the  considered  oceanic  data  were  acquired 
between  1981  and  2006.  Most  importantly,  this  period 
encompasses  all  the  surface  drifter  observations  which  have  a 
major  impact  on  the  reconstructed  MDT,  and  coincides  with 
the  time  span  of  massive  satellite  observations  of  the  ocean 
surface.  For  these  reasons  we  assume  that  the  reconstructed 
MDT  represents  best  the  20  year  period  of  1986-2006. 

2,2.  Analysis  Technique 

[  19]  The  mean  climatological  state  of  the  Bering  Sea  was 
found  as  a  data-optimized  solution  of  the  primitive  equation 
model  already  used  by  Panteleev  el  al.  [2006b]  for 
retrieving  the  mean  summer  circulation  in  the  northern  part 
of  the  basin  and  in  the  Kara  Sea  [Panteleev  el  al .,  2007].  The 
numerical  model  is  a  modification  of  the  C  grid,  z  coordinate 
Ocean  General  Circulation  Model  (OGCM)  designed  by 


Madec  et  al  [1999]  (see  Nechaev  el  al.  [2005]  for  details). 
The  model  was  configured  in  the  domain  shown  in  Figure  1 
with  meridional  and  zonal  resolutions  of  0.16°  and  0.3°, 
respectively.  The  corresponding  grid  step  (18  km)  is  larger 
than  the  Rossby  deformation  radius  to  suppress  mcsoscale 
eddies,  but  it  is  fine  enough  to  resolve  the  major  topographic 
and  circulation  features,  including  the  Near  Strait,  the 
Amchitka  and  Amukta  passes,  and  the  Kamchatka  Current. 
Vertically,  the  grid  had  34  levels  with  spacing  ranging  from 
5  m  near  the  surface  to  500  m  in  the  deep  layers. 

[20]  The  model  resolution  (18  km)  was  chosen  to  be 
somewhat  larger  than  typical  resolution  of  the  altimetry 
observations  (10  km),  which  barely  resolve  the  local  defor¬ 
mation  radius  (15  km)  and  mesoscale  eddy  dynamics.  Apart 
from  the  immense  computational  cost  of  the  eddy-resolving 
4DVar  assimilation  [c.g.,  Ho  tell  et  al ,  2010;  Mazloff  et  al, 
2010],  this  choice  of  a  coarser  grid  was  made  to  avoid  reg¬ 
ularization  of  the  adjoint  model  through  the  artificial  increase 
of  the  horizontal  diffusion.  Effectively,  such  an  approach 
performs  optimization  in  the  subspacc  of  smooth  model 
solutions,  consistent  with  the  large  diffusion  of  the  adjoint 
model,  and  may  not  converge  on  the  “true"  optimal  solution 
[c.g.,  Yaremchuk  et  al,  2009]. 

[21]  A  4DVar  data  assimilation  algorithm  was  configured 
to  find  a  quasi -stationary  solution  to  model  equations  that 
optimally  fits  the  oceanic  and  atmospherical  observations 
described  above.  Following  the  approach  of  Tziperman  and 
Thacker  [1989],  in  the  first  guess  we  specified  initial  and 
boundary'  conditions  and  integrated  the  model  for  a  period 
of  several  weeks  w  ith  steady  climatological  momentum  and 
heat/salt  fluxes  at  the  surface.  Both  surface  forcing  and 
initial/boundary  conditions  were  imposed  as  weak  con¬ 
straints,  i.e.,  they  were  iteratively  optimized  in  the  course  of 
assimilation.  Technically,  the  optimization  procedure  can  be 
viewed  as  a  dynamically  constrained  minimization  of  the 
cost  function  J  which  measures  the  distance  between  the 
model  solution  and  the  data.  The  minimization  is  performed 
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Figure  2.  (a)  Preproeessed  temperature  data  at  0  m  used  in  assimilation  and  (b)  the  associated  standard 
deviations  at  0  m.  Contour  arc  in  °C. 


using  the  standard  Lagrangian  multiplier  technique  [e.g., 
Le  Dimet  and  Talagrand ,  1986]  by  adjusting  the  errors  in 
the  fields  that  directly  control  the  model  solution.  These  error 
fields  eF  included  differences  between  atmospheric  forcing 
fields  F  and  their  available  estimates  F*,  as  well  as  similar 
differences  between  the  initial/open  boundary  conditions  for 
temperature,  salinity,  horizontal  velocity  v  =  {w,  v}  and  sea 
surface  height  £.  The  total  number  of  the  grid  points  occu¬ 
pied  by  the  control  fields  (the  dimension  of  the  control 
vector  c)  was  N  =  856,054.  In  addition,  the  cost  function 
included  smoothness  constraints,  which  penalized  grid-scale 
components  of  the  model  fields,  and  the  “steady  state  con¬ 
straint,”  which  enforced  quasi-stationary  state  by  penalizing 
the  difference  between  the  model  fields  at  the  beginning  and 
at  the  end  of  integration. 


[22]  The  cost  function  has  the  form 


r  [Y[0)  -  Y(T)]2— 


+  j  »'f(AK)2^  +  J eFf)j. 

n  s 


(1) 


Here  Y  -  {T,  S,  v,  denotes  the  set  of  prognostic  model 
fields  and  corresponding  data,  F  stands  for  the  set  of  surface 
forcing  fields,  A  is  the  Laplacian  operator,  Ft  is  the  3D 
model  domain,  S  is  its  upper  boundary  at  2  =  0,  and  the 
overbar  denotes  time  average  over  the  period  of  integration. 
Summation  in  the  first  term  is  made  over  the  P  =  403,156 
observation  points  described  in  section  2.1.  The  index  p 
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Figure  3.  (a)  The  multiyear  mean  surface  drifter  velocities  averaged  over  the  model  grid  cells.  The  gray 
scale  shows  the  RMS  variance  (in  cm/s)  of  the  drifter  velocity  in  each  grid  cell,  (b)  Long-term  mooring 
velocities  (gray  arrows)  in  the  upper  200  m  used  in  assimilation  and  averaged  velocities  of  Argo  drifters  at 
1000  m  (not  used  in  assimilation).  Several  typical  trajectories  of  the  Argo  drifters  arc  shown. 


enumerates  spatial  locations  of  the  observations  Y*  and  the 
components  of  Y  observed  at  these  locations. 

[23]  Weighting  functions  W  before  the  squared  quantities 
have  the  sense  of  the  inverse  error  variances,  so  that  the  cost 
function  can  be  interpreted  as  an  argument  of  the  multi¬ 
variate  Gaussian  probability  distribution  [ Thacker ,  1989]. 
Linder  such  interpretation  the  optimal  solution  is  the  most 
probable  model  state  for  a  given  set  of  observations  over  a 
specific  time  period  and  their  prior  error  statistics. 

2.3.  Error  Estimation 

[24]  Since  the  cost  function  (1)  implicitly  depends  on  the 
set  of  the  adjusted  parameters  c  (control  variables),  a  pos¬ 


teriori  error  eo variance  can  in  principle  be  obtained  as  the 
inverse  of  the  Hessian  matrix  H  =  c^J/dc2  [Thacker,  1989]. 

[25]  In  practice,  inversion  of  the  N  *  N  Hessian  matrix  is 
computationally  prohibitive,  so  we  employed  an  approxi¬ 
mate  approach  assuming  that  the  structure  of  the  SSH  error 
fields  roughly  followed  their  patterns  of  natural  variability. 
These  patterns  were  estimated  as  the  leading  empirical 
orthogonal  functions  (EOFs)  of  the  covariance  matrix  C 
derived  from  the  Aviso  SSH  anomalies.  The  corresponding 
spectrum  shows  that  90%  of  the  SSH  variability  (and  error 
variance)  could  be  explained  by  n  =*  20  modes  (Figure  4). 
The  Hessian  matrix  was  inverted  in  the  n -dimensional 
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Figure  4.  Normalized  spectrum  of  the  SSH  covariance 
matrix  (solid  curve)  derived  from  the  Aviso  SSH  anomalies 
regridded  on  the  model  domain.  The  dashed  curve  shows 
percentage  of  the  SSH  variance  unexplained  by  the  given 
number  of  modes  n. 

subspace  spanned  by  these  leading  modes  and  then  pro¬ 
jected  back  on  the  model  grid  using  the  relationship 

Cc  =  P(P7HP)~'pr.  (2) 

Here  P  is  the  20  *  8,722  matrix  whose  columns  are  the  20 
eigenvectors  of  C  with  the  largest  eigenvalues,  8,722  is  the 
number  of  SSH  model  field  grid  points,  and  T  denotes 
transposition.  Mapping  the  square  root  of  the  diagonal  ele¬ 
ments  of  C^-  gives  an  estimate  of  the  uncertainty  of  the 
optimized  MDT  (see  section  3.1). 

[2f>]  The  error  covariance  of  the  surface  gcostrophic 
velocities  is  estimated  using  the  relationship 

Cv  =  GyrCcGv,  (3) 

where 


is  the  matrix  composed  of  the  finite  difference  representa¬ 
tions  of  the  operators  dY  and  ~dx.  Here  g  is  the  gravity 
acceleration,  -  1 .025  g/cm1  is  the  mean  density  of  sea¬ 
water  and/is  the  Coriolis  parameter.  The  spatial  distribution 
of  the  velocity  error  variance  is  obtained  by  taking  the 
square  root  of  the  diagonal  elements  of  Cv. 

[27]  Seasonality  of  the  T/S  observations  may  result  in 
biasing  the  MDT  toward  summer  conditions.  To  estimate 
this  effect  we  calculated  separately  the  baroclinic  impacts  of 
the  winter  and  summer  T/S  distributions  on  the  MDT.  It  was 
found  that  the  mean  RMS  difference  between  such  “sea¬ 
sonal”  MDT  estimates  is  about  3  cm.  The  respective  MDT 
distributions  show  differences  at  small  spatial  scales  in  the 
major  part  of  the  domain.  The  only  exception  is  the  Alaskan 
shelf  in  the  east  where  the  difference  between  the  summer 


and  winter  MDTs  is  the  most  profound.  Overall,  wc  estimate 
that  seasonality  in  the  T/S  data  may  result  in  a  bias  of  1  -2  cm 
over  the  major  part  of  the  basin  with  an  increase  to  2-3  cm  on 
the  Alaskan  shelf. 

3.  MDT  of  the  Bering  Sea 

3.1.  General  Features 

[2s]  Since  the  optimized  surface  velocities  below  the 
Ekman  layer  are  in  approximate  geostrophic  balance  with 
the  MDT  (optimized  SSH  field),  the  reconstructed  MDT 
contours  (Figure  5a)  are  conveniently  interpreted  as 
streamlines  of  the  mean  geostrophic  currents  at  the  surface. 
The  circulation  pattern  reveals  the  following  major  struc¬ 
tures:  (1)  an  intense  (30-40  cm/s)  Alaskan  Stream  south  of 
the  Alaska  Peninsula,  (2)  a  somewhat  weaker  (10-20  cm/s) 
Aleutian  North  Slope  Current  embracing  the  southern  and 
northern  flanks  of  the  Aleutian  Arc,  (3)  the  30-40  cm/s 
strong  Kamchatka  Current  on  the  west,  and  (4)  a  relatively 
weak  (5-1 5  cm/s)  cyclonic  circulation  occupying  the  deep 
part  of  the  Bering  Sea.  According  to  Figure  5a,  a  significant 
portion  of  this  cyclonic  gyre  originates  in  the  Near  Strait 
with  the  rest  coming  from  the  inflow  through  other  Aleutian 
passages. 

[29]  The  circulation  shown  in  Figure  5a  is  in  good  qual¬ 
itative  agreement  with  the  results  of  Stabeno  et  al.  [2005] 
who  describe  gradual  leakage  of  the  Alaskan  Stream 
though  the  passages  in  the  Aleutian  Arc.  In  their  earlier 
work,  Stabeno  and  Reed  [1994]  estimated  the  splitting  point 
of  the  Bering  Slope  Current,  which  occupies  the  eastern 
flank  of  the  deep  cyclonic  gyre.  According  to  their  analysis 
of  surface  drifters,  the  current  splits  into  two  branches 
around  60°N,  1 76°E:  at  this  point  a  larger  part  of  the  surface 
flow  begins  to  form  the  Kamchatka  Current,  while  the  rest 
takes  the  path  to  the  Arctic  Ocean.  In  Figure  5a  this  splitting 
occurs  at  61°N  1 78°W,  an  insignificant  difference  given  the 
2-3  cm  uncertainty  of  our  result  in  this  area  (Figure  6). 

[30]  The  pattern  in  Figure  5a  docs  not  reveal  a  clear  split 
of  the  Kamchatka  Current  into  coastal  and  offshore  branches 
in  the  vicinity  of  the  Shirshov  Ridge  as  was  shown  by 
Panteleev  et  al.  [2006b].  Instead,  a  relatively  broad  west¬ 
ward  current  across  the  Kamchatka  Basin  has  been  obtained. 
This  difference  can  be  attributed  to  the  seasonal  nature  of 
the  “offshore  branch”  of  the  Kamchatka  Current. 

[3 1  ]  The  mean  mismatch  betw  ecn  the  reconstructed  surface 
velocities  and  assimilated  drifter  velocities  (Figure  3a)  is 
8.2  cm/s.  It  is  unlikely  that  better  agreement  could  be 
obtained  between  the  climatological  velocity  with  mean 
amplitude  of  approximately  10  cm/s  and  the  highly  variable 
surface  currents  derived  from  the  drifter  trajectories  affected 
by  eddies  and  small-scale  variations  of  the  wind  stress. 
Similar  discrepancies  (on  the  order  of  7-9  cm/s)  were 
observed  by  Maximenko  et  al  [2009]  who  estimated  velocity 
errors  in  three  global  MDT  products  and  by  Panteleev  et  al 
[2006b]  who  also  obtained  similarly  high  model-surface 
drifter  velocity  misfits.  Despite  this,  the  optimized  surface 
velocities  in  the  core  of  the  Kamchatka  Current  are  about 
22  cm/s,  in  relatively  good  agreement  with  local  drifter 
velocities  (20-30  cm/s)  [sec  Hughes  et  al ,  1974]  (also  see 
Figure  3a).  Similar  agreement  was  obtained  for  the  Alaskan 
Stream.  For  example,  Stabeno  et  al  [2005]  observed  a 
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Figure  5.  The  mean  dynamic  topographies  of  the  Bering  Sea  obtained  (a)  in  the  present  study,  (b)  by  Rio 
et  al.  [2005],  (c)  by  Rio  et  al.  [2009],  and  (d)  by  merging  the  EGM08  geoid  model  with  altimeter  data. 
Validation  subdomains  for  the  presented  MDT  are  shown  by  solid  black  rectangles  for  drifter  inter¬ 
comparison  experiments  and  by  white  rectangle  for  4DVar  intercomparison  experiment.  Contour  interval 
is  4  cm. 


E  165  170  175  180  -175  -170  -165  -160  W 

Figure  6.  A  posteriori  RMS  error  variance  map  for  the  MDT.  Contour  interval  is  2  cm. 
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4  year  mean  velocity  of  40  cm/s  at  100  m  near  the  Amukta 
Pass,  which  is  close  to  our  optimized  estimate  of  30  cm/s. 

[3:]  To  further  validate  the  MDT,  we  compared  the 
optimized  velocity  field  at  1000  m  with  independent 
velocity  data  derived  from  the  Argo  floats.  At  this  depth  the 
average  velocity  of  the  floats  is  4.6  cm/s,  which  is  close  to 
the  mean  optimized  velocity  magnitude  of  3.7  cm/s.  We 
consider  this  to  be  in  good  agreement  with  observations 
because  a  significant  fraction  of  Argo  drifters  are  involved 
in  mesoscale  motions  resulting  in  a  higher  mean  Lagrangian 
velocity  compared  to  the  mean  Eulerian  velocity.  Inter¬ 
estingly,  the  optimized  velocities  in  the  Kamchatka  Current 
(10  cm/s)  and  the  Alaskan  Stream  (8  cm/s)  at  1000  m 
agree  almost  perfectly  with  the  estimates  of  1 1  cm/s  and 
8  cm/s  derived  from  Argo  drifters  (Figure  3b). 

[33]  In  the  major  Aleutian  Passes  the  volume  transports 
driven  by  the  MDT  shown  in  Figure  5a  were  found  to  be 
2.5-7  times  larger  than  the  those  obtained  by  Stabeno  ct  al 
[1999]  by  the  dynamical  method.  This  discrepancy  is  likely 
due  to  underestimation  of  the  barotropic  velocity  in  the 
Aleutian  straits  by  the  dynamical  method.  Recent  observa¬ 
tions  by  Stabeno  ct  al.  [2005]  seem  to  support  this  sug¬ 
gestion,  as  they  revealed  substantial  (up  to  50  cm/s) 
northward  currents  at  100-200  m,  indicating  similarity  with 
our  result  and  the  barotropic  nature  of  the  flow  through  the 
Aleutian  passes. 

[34]  Inspection  of  the  error  map  in  Figure  6  shows  that  the 
MDT  errors  are  significantly  smaller  (2-3  cm)  than  the 
typical  error  estimates  of  6-8  cm  obtained  for  the  global 
products  [e.g.,  Maximenko  et  al.,  2009].  Such  improvement 
can  be  explained  by  the  extensive  database  and  stronger 
dynamical  constraints  underlying  the  new  MDT.  The  pro¬ 
nounced  maximum  in  the  error  field  between  62°  and  64°N 
(Figure  6)  is  caused  by  ice  cover  obstructing  observations  in 
winter.  A  certain  decrease  of  the  uncertainty  north  of  St. 
Lawrence  Island  is  due  to  the  relatively  dense  mooring 
observations  in  that  area  (cf.  Figure  3b).  The  respective 
velocity  error  v  ariances,  computed  using  (3)— (4),  range  from 
2-4  cm/s  indicating  robustness  of  the  flow  pattern  shown  in 
Figure  5a. 

[35]  It  is  necessary  to  note  that  our  approach  does  not  take 
into  account  residual  tidal  velocities,  which  do  not  exceed 
0.2-0.5  cm/s  in  the  open  sea,  but  can  be  as  large  as  2-3  cm/s 
along  the  continental  slope  and  reach  5-10  cm/s  around  the 
islands  of  the  Aleutian  Arc  [Kowalik,  1999].  At  the  same 
time,  the  residual  tidal  velocities  arc  known  to  create  trapped 
clockwise  circulations  around  the  islands,  thus  minimizing 
their  effect  on  the  total  volume  transport  in  the  open  Bering 
Sea. 

[3ft]  Overall,  the  obtained  MDT  and  associated  currents 
are  in  good  agreement  with  the  previous  studies  and  allow 
reasonable  quantitative  assessment  of  the  mean  background 
circulation  in  the  Bering  Sea.  Importantly,  the  SSH  pattern 
shown  in  Figure  5a  is  dynamically  balanced  with  the  cli¬ 
matological  temperature  and  salinity  fields  in  the  region.  In 
section  3.2  we  validate  the  product  against  MDTs  obtained 
by  alternative  methods. 

3.2.  Validation  Against  Other  Products 
3.2.1.  MDTs  of  the  Bering  Sea 

[37]  An  unprecedented  increase  in  the  amount  of  data  on 
surface  winds,  currents,  and  SSH  anomalies  has  fueled 


numerous  efforts  to  obtain  accurate  MDT  estimates  from 
observations  constrained  by  simplified  dynamics.  Niileret  al. 
[2003]  employed  surface  drifter  velocities,  wind  stress,  and 
SSH  anomalies  to  produce  the  first  global  MDT  at  0.5° 
resolution.  In  a  more  comprehensive  effort,  Rio  et  al.  [2005, 
hereinafter  R05]  combined  these  data  with  temperature  and 
salinity  from  the  World  Ocean  Atlas  [Conkright  et  al. ,  2002] 
to  produce  a  more  accurate  estimate.  Most  recent  efforts  by 
Maximenko  et  al.  [2009]  and  Rio  et  al.  [2009,  hereinafter 
R09]  employ  diverse  in  situ  and  remotely  sensed  observa¬ 
tions  in  conjunction  with  the  GGM02C  gravity  model 
[Riegber  et  al.,  2005].  A  good  MDT  estimate  based  on 
the  Earth  Gravitational  Model  2008  (EGM08)  geoid  and 
altimeter  data  only  was  obtained  recently  by  Andersen 
and  Knudsen  [2009,  hereinafter  A09]  (see  also  http:// 
gracetellus.jpl.nasa.gov/data/dot/). 

[38]  All  these  efforts  were  made  on  the  global  scale  with 
resolutions  0.5°  or  higher.  Additional  difficulties  in 
retrieving  the  MDT  from  observations  do  arise  in  subpolar 
regions,  where  drifters  and  altimetry  are  obscured  by  icc  and 
the  dominant  scales  get  smaller  due  to  a  decreasing  defor¬ 
mation  radius. 

[39]  Figure  5  gives  a  comparison  of  our  product  with  three 
latest  MDTs:  two  developed  by  Rio  and  others  in  the  Aviso 
altimetry  processing  center  (Figures  5b  and  5c)  and  the  third 
one  by  the  Danish  Space  Center  and  the  Jet  Propulsion 
Laboratory'  (Figure  5d).  All  four  patterns  have  much  in 
common:  the  Alaskan  Stream  south  of  Alaska  Peninsula, 
northeasterly  Slope  Current  in  the  center,  and  the  Kam¬ 
chatka  Current  on  the  east.  Figures  5b-5d  also  demonstrate 
a  basin-scale  SSH  difference  of  40-50  cm  between  the 
Alaska  and  the  southwestern  comer  of  the  domain,  which  is 
15  cm  less  than  in  Figure  5a.  There  arc,  however,  much 
larger  quantitative  differences  between  the  patterns.  First, 
the  Kamchatka  Current  is  barely  visible  in  Figure  5c,  and 
especially  in  Figures  5b  and  5d  where  SSH  contours  have 
the  correct  offshore  slope  only  along  limited  portions  of 
the  coastline.  R09  is  more  realistic,  but  they  still  under¬ 
estimate  the  current’s  width  and  velocity.  Second,  MDTs 
in  Figures  5b-5d  drop  12-18  cm  between  St.  Lawrence 
island  and  the  Bering  Strait,  correctly  indicating  that  the 
Bering  Strait  transport  is  driven  by  the  SSH  difference 
between  the  Pacific  and  the  Arctic  oceans,  but  the  flow 
through  the  strait  appears  to  be  completely  agcostrophic.  In 
contrast,  the  presented  MDT  finds  a  25  em  difference 
between  St.  Lawrence  Island  and  the  Chukchi  Peninsula  and 
the  mean  flow  through  the  Bering  Strait  is  geostrophically 
balanced.  This  agrees  well  with  the  results  of  Cherniawsk)' 
et  al.  [2005],  who  found  that  the  Bering  Strait  transport 
can  be  retrieved  from  sea  level  anomalies  with  reasonable 
accuracy  using  gcostrophy.  Finally,  the  cyclonic  gyre  in  the 
deeper  southwestern  part  of  the  Bering  Sea  is  not  clearly 
visible  in  the  R05  and  R09  MDT  maps.  Existence  of  this 
gyre  is  supported  by  drifter  data  (Figure  3a)  [Stabeno  et  al., 
1999;  Johnson  et  al.,  2004],  data-constrained  model  simu¬ 
lations  [Yaremchuk,  2001;  Wunsch  et  al,  2009];  and  the 
A09  MDT  (Figure  5d)  which  is  presumably  more  accurate 
in  the  open  sea.  The  A09  pattern  also  appears  to  be  more 
consistent  with  the  Argo  float  data  (Figure  3b)  and  the  result 
of  Stabeno  and  Reed  [1994]  as  it  shows  a  clear  meridional 
inflow  into  the  domain  at  1 67°  1 70°E. 
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Table  1.  Difference  (cm/s)  Between  the  Geostrophic  Currents 
Obtained  With  Various  MDT  Products  and  Drifter  Velocity  Data 
Used  in  the  Present  Study1* 


Region 

R05 

A09 

R09 

4DVar 

KC 

12.0 

13.2 

11.0 

9.7 

BSC 

10.1 

9.5 

9.6 

9.1 

AS 

17.6 

17.3 

18.0 

15.4 

Mean 

13.2 

13.3 

12.8 

11.4 

"Most  of  lhc  drifters  were  drogued  at  40  m 


[40]  Qualitative  comparison  of  Figures  5a-5d  emphasizes 
the  validity  of  additional  in  situ  data  in  reconstructing  the 
MDT  in  nearshore  regions,  where  altimetry  and  gravity 
observations  tend  to  be  less  aeeurate  (ef.  Figure  5d  and 
Figures  5a-5e  near  the  Aleutian  are).  On  the  other  hand,  the 
EGM08  geoid  appears  to  be  more  aeeurate  in  the  eenter  of 
the  Sea,  as  it  is  able  to  reproduce  the  eyelonie  gyre  over  the 
deep  part  of  the  basin  and  the  above  mentioned  meridional 
inflow'  around  1 70°E. 

3.2.2.  Quantitative  V  alidation 

[41]  To  assess  the  MDT  products  quantitatively,  we  used 
the  Aviso  methodology:  three  domains  well  eovered  by 
drifters  were  pieked,  gridded  Aviso  SSH  anomalies  from 
http://www.aviso.oeeanobs.eom  were  added  to  eaeh  MDT, 
and  the  resulting  geostrophic  currents  were  compared  with 
the  currents  dedueed  from  drifter  trajectories.  The  test 
domains  (Figure  5a)  eover  the  major  circulation  features  of 
the  Bering  Sea  and  include  the  Kamchatka  Current  (KC), 
the  Bering  Slope  Current  (BSC),  and  the  Alaskan  Stream 
(AS).  Table  1  summarizes  the  calculated  RMS  velocity 
discrepancies. 

[42]  The  presented  4DVar  MDT  demonstrates  smaller 
errors,  especially  in  the  KC  region,  where  all  the  alternative 
MDTs  fail  to  reproduce  a  continuous  eurrent  with  realistic 
velocities.  This  sueeess  may  be  partially  due  to  the  faet  that 
the  mean  drifter  velocities  (Figure  3a)  were  assimilated  in 
the  4DVar  MDT.  Note  that  both  the  R05  and  the  R09 
produets  were  also  developed  by  taking  the  drifter  infor¬ 
mation  into  aeeount.  Table  1  also  shows  gradual  improve¬ 


ment  of  the  MDT  products  with  time:  there  is  an  overall 
decrease  of  the  mean  error  from  left  to  right  indieating  a 
tendency  of  convergence  to  the  “true  MDT.“  It  is  also 
noteworthy  that  A09  outperforms  R09  in  the  BSC  and  AS 
regions;  although  it  was  compiled  with  altimetry  data  only, 
it  utilized  a  better  geoid. 

[43]  In  addition,  we  validated  the  performance  of  the 
above  MDTs  in  approximating  the  nearshore  geostrophie 
currents  in  Kuskokwim  Bay,  where  32  drifters  were  laun¬ 
ched  by  Quinhagak  fishermen  in  July-Oetober  2008. 

[44]  Using  the  4DVar  algorithm  described  in  seetion  2, 
drifter  velocities  were  assimilated  together  with  elimato- 
logieal  T/S  distributions  into  the  model  configured  in  the 
40  x  40  domain  shown  in  Figure  7.  The  optimized  SSH  field 
was  averaged  over  the  assimilation  period  and  used  as  a 
benchmark  for  assessing  the  quality  of  the  MDTs.  This 
assessment  was  done  by  comparing  the  above  mentioned 
time-averaged  distribution  Q,  (Figure  7a)  with  the  sum  of 
MDT  and  Aviso  anomalies  averaged  over  the  period  of  the 
Kuskokwim  experiment. 

[45]  To  quantify  the  assessment,  we  also  computed  rela¬ 
tive  errors 

£=\/<(0-002>/<( 

where  angular  brackets  stand  for  the  average  over  the  region 
with  drifter  trajectories  shown  by  the  white  rectangle  in 
Figure  7a.  The  values  of  were  found  to  be  0.79,  1.71, 
1.06,  and  1.27  for  the  4DVar,  R05,  R09,  and  A09  MDTs, 
respectively.  These  large  values  are  explained  by  higher 
MDT  errors  in  shallow'  regions,  errors  in  Aviso  SSH 
anomalies,  and  errors  in  Q>  itself.  Nevertheless,  the  mean 
SSH  distribution,  defined  as  the  sum  of  the  Aviso  anomalies 
and  4DVar  MDT  (Figure  7b),  is  in  better  agreement  (£  = 
0.79)  with  the  benchmark  field  (Figure  7a)  than  with  the 
corresponding  SSH  distributions  derived  from  the  R05 
(Figure  7c,  £  =  1.71),  A09  (£  =  1.27),  or  R09  (£  =  1.09) 
MDTs. 

[46]  All  these  examples  illustrate  the  significant  advantage 
of  the  presented  MDT  as  compared  to  the  best  global  pro- 


Figure  7.  (a)  Mean  SSH  derived  from  assimilating  drifter  velocities  in  the  Kuskokwim  Bay  and  the 
mean  SSH  obtained  by  adding  the  (b)  presented  MDT  and  (e)  MDT  of  Rio  et  al.  [2005]  to  Aviso 
SSH  anomalies  averaged  over  the  duration  of  the  Kuskokwim  drifter  experiment  (7  July  2008  to  12 
October  2008). 
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ducts,  which  generally  suffer  from  larger  errors  in  near- 
coastal  regions.  Apart  from  somewhat  better  resolution,  our 
product  benefits  from  much  better  usage  of  the  dynamical 
information  which  significantly  constrains  the  SSH  field 
near  the  coast. 

4.  Summary  and  Discussion 

[47]  A  quantitative  estimate  of  the  MDT  in  the  Bering  Sea 
has  been  presented.  The  product  is  obtained  by  combining 
historical  oceanographic  and  atmospheric  observations  with 
high-resolution  model  dynamics  in  the  framework  of  the 
variational  lechnique.  The  optimized  circulation  is  dynami¬ 
cally  balanced  and  statistically  consistent  with  the  utilized 
data. 

[48]  The  presented  dynamic  topography  is  validated 
against  the  three  most  recent  global  MDT  products.  Two 
of  them  (R05  and  R09)  were  derived  by  merging  in  situ 
observations  with  winds,  satellite  altimetry,  and  gcoid 
models.  The  third  product  (A09)  is  based  only  on  altimeter 
data  and  the  most  recent  EGM08  geoid.  Comparison  has 
shown  that  the  first  two  MDTs  (R05  and  R09)  perform 
better  than  A09  near  the  land,  but  are  nevertheless  notice¬ 
ably  worse  than  the  presented  product.  On  the  other  hand, 
A09  docs  much  better  job  in  representing  the  SSH  pattern  in 
the  open  Bering  Sea,  indicating  a  considerable  increase  in 
the  accuracy  of  the  EGM08  geoid  compared  to  the  previous 
models.  Overall,  the  presented  MDT  has  shown  better  per¬ 
formance  against  independent  drifter  and  mooring  observa¬ 
tions  while  demonstrating  somewhat  more  realistic 
geostrophic  circulation,  especially  near  the  Kamchatka 
coastline  and  north  of  St.  Lawrence  island. 

[49]  We  conducted  a  rigorous  error  analysis  of  the  MDT 
and  related  geostrophic  currents.  The  MDT  uncertainties 
were  computed  by  inverting  the  Hessian  matrix  of  the 
assimilation  problem  in  the  subspacc  spanned  by  the  gravest 
modes  of  SSH  variability  observed  from  satellites.  Geo¬ 
strophic  velocity  errors  were  calculated  by  the  appropriate 
transformation  of  the  SSH  error  covariance  matrix.  This 
analysis  has  shown  the  remarkable  robustness  of  our  MDT 
estimate,  which  is  characterized  by  SSH  errors  of  2-3  cm, 
considerably  less  than  typical  error  variances  (6-8  cm)  of 
the  global  MDTs. 

[50]  The  corresponding  geostrophic  velocity  errors  range 
within  2-4  cm/s  making  it  possible  to  quantify  the  major 
circulation  features  in  the  Bering  Sea  with  a  reasonable 
degree  of  confidence.  In  particular,  the  mean  surface  cur¬ 
rents  in  the  AS  (25  40  cm/s)  and  KC  (15-30  cm/s)  arc 
determined  with  a  formal  accuracy  of  10-15%,  and  in  the 
Aleutian  North  Slope  Current  ( 1 0—20  em/s)  with  an  accu¬ 
racy  of  20-25%. 

[51]  Our  comparisons  with  similar  products  also  indicate 
that  the  A09  MDT  estimate  derived  from  the  latest  GRACE 
gcoid  (EGM08)  has  more  in  common  with  the  presented  sea 
surface  topography  than  the  MDTs  based  on  earlier  versions 
of  the  gcoid,  especially  in  the  open  sea  regions,  where 
altimeter  observations  tend  to  be  more  accurate. 

[52]  Having  a  realistic  SSH  reference  is  especially 
important  for  successful  monitoring  of  the  Bering  Sea  cir¬ 
culation.  Currently,  the  Jason,  Envisat,  and  GFO  satellite 
altimeter  missions  provide  accurate  SSH  anomalies  across 
the  entire  Bering  Sea  every  10-30  days,  but  because  of  the 


insufficient  knowledge  of  MDT/gcoid,  use  of  these  data  is 
not  straightforward.  We  believe  that  the  presented  MDT  will 
improve  the  accuracy  in  estimating  the  Bering  Sea  surface 
circulation  and  may  even  be  used  for  calibrating  the  gcoid 
models  in  the  region. 

[53]  Recent  changes  in  the  Arctic  and  North  Pacific  Cli¬ 
mate  [e.g.,  Weller ,  1998;  Oxerland  et  al.,  1999,  2000]  may 
seem  to  contradict  the  validity  of  the  steady  state  assumption 
underlying  our  analysis.  This  assumption,  however,  is 
imposed  in  the  form  of  a  weak  constraint  (second  term  in 
equation  (1))  and  allows  residual  trends  in  the  T/S  fields. 
Their  magnitudes  of  0.15°C/yr  appear  to  be  consistent  with 
the  existing  experimental  estimates  for  the  period  of  1 974— 
1994  described  by  Luchin  et  al  [2002]. 

[54]  The  presented  MDT  can  be  viewed  and  downloaded 
from  http://people.iarc.uaf.edu/^gleb/nprb_aleutian_passes/ 
bering  sea  atlas  rcgister.php  together  with  the  optimized 
climatological  temperature,  salinity,  and  velocity  distribu¬ 
tions.  These  fields  were  obtained  simultaneously  with  the 
4DVar  reconstruction  of  the  MDT  but  error  variance  dis¬ 
tributions  have  not  yet  been  produced  for  them.  Wc  are 
currently  working  on  improving  error  estimates  for  the  MDT 
and  the  other  fields  of  this  climatological  atlas. 

[55]  The  proposed  approach  is  a  relatively  inexpensive 
way  to  use  diverse  observational  data  in  deriving  an  MDT 
for  any  region.  Because  the  period  of  model  integration  is 
relatively  short,  the  method  can  also  be  viewed  as  an 
“iterative  diagnostic  calculation"  with  updated  initial  and 
boundary  conditions.  Wc  compared  the  4DVar  MDT  in  the 
southeastern  part  of  our  domain  with  the  MDT  reccnlly 
proposed  for  the  Gulf  of  Alaska  by  Foreman  et  al.  [2008], 
and  found  a  very  good  agreement  between  these  products 
for  the  AS  region.  That  shows  that  these  two  dynamically 
constrained  approaches  have  much  in  common.  The  limiting 
factor  for  our  approach  is  model  resolution,  w  hich  should  be 
kept  relatively  low,  while  the  diagnostic  calculation  can  be 
formally  done  at  a  very  fine  resolution  of  1  km  [Foreman  et 
al .,  2008].  The  variational  approach  allows  us,  however,  to 
take  into  account  all  available  data  together  with  ihcir  sta¬ 
tistics  and  provides  a  consistent  formalism  for  estimating  a 
posteriori  errors. 
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